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We present integral field spectroscopic observations of the central region of the 
active galaxy NGC 4258 obtained with the fibre IFU system INTEGRAL. We have 
been able to detect cold neutral gas by means of the interstellar NaD doublet absorp- 
tion and to trace its distribution and kinematics with respect to the underlying disc. 
The neutral gas is blue-shifted with projected velocities in the 120-370 km/s range. 
We have also detected peculiar kinematics in part of the ionized gas in this region 
by means of a careful kinematic decomposition. The bipolar spatial distribution of 
the broader component is roughly coincident with the morphology of the X-ray dif- 
fuse emission. The kinematics of this gas can be explained in terms of expansion at 
very high (projected) velocities of up to 300 km/s. The observations also reveal the 
existence of a strip of neutral gas, parallel to the major kinematic axis, that is nearly 
coincident with a region of very high [SII]/Ha ratio tracing the shocked gas. Our 



observations are consistent with the jet model presented by IWilson et al.l (|2001f ) in 
which a cocoon originating from the nuclear jet is shocking the gas in the galaxy disc. 
A lternat i vely, our o bservations are also co nsistent with the bipolar hypershell model 
of lSofud (|1980l ) and lSofue fc Voglerl (|200hY On balance, we prefer the latter model as 
the most likely explanation for the puzzling features of this peculiar object. 
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1 INTRODUCTION 

The nearby spiral galaxy NGC 4258 (M106) is a particu- 
larly interesting galaxy. It presents two peculiarities that 
have been the cause of many studies of this object: the con- 
firmation of the existence at its centre of a massive black 
hole, detected in 1995 through measurements o f high speed 
motions in the water masers (jMivoshi et al.1 |l995)), and 
the existence of t he mysterious anomalous ar ms discovered 
in the optical by [Courtes fc Cruvellierl (| 196 lT ) and in radio 
by Ivan der Kruit et al.l l|l972[ ). As we shall see below, it is 
indeed not unlikely that both phenomena are closely related. 

These anomalous arms are two large arm-like gaseous 
structures of uncertain origin that extend through most 
of the visible disc of the galaxy. There have been several 
possible explanations for the origin of these arms, most 
of which can be traced back to three distinct possibili- 
ties: (a) expelled gas from the nucleus in t h e plane, ei- 
t her i n a ballistic Ivan der Kruit et al.l d 19721) . Ivan Albadal 
1978), 



van Albada fc van der Hu lst (1982)) or a jet model 



.,Ford et al.1 dl98rj), ICecil et al.1 (119921 )), (b) a bar-induced 
shock front (|Cox fc Downed ( 19961 )). and (c) gas expelled out 



of the p l ane ei t her in a jet, hypershell , etc. (|Sofuel l|l980j ). 
ISandersI l|l982l ), ISofue fc Voglerl <|200ll )). There are indeed 
many observational features that have to be explained by 
an y suitable model (see for example the exc ellent works 
by ICecil et al.l (|2000l ) and IWilson et all (|200ll ) including a 
plethora of optical, radio and X-ray data) . 

The morphology of the anomalous arms can be greatly 
influenced by projection effects. Although it is pretty well 
established that the inner parts of these arms are in the disc 
plane, it is not known for certain whether the arms are three 
dimensional structures seen in projection (and enhanced by 
a limb-brightening effect), or if they are indeed mostly con- 
tained in a plane (and if this is the case, the orientation of 
this plane). On the other hand, the highly distorted kine- 
matics that appear in the innermost region of the galaxy 
also imposes restrictions on the possible scenarios. 

The galaxy NGC 4258 has an inclination of about 60- 
64°, and a position ang le of 146-160° (cf. IJarrett et al.l 
(|2003t) . ICecil et all l|l992l )). The strong absorption below the 
major axis shows clearly that the region located SW of the 
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major axis is the near side of the galaxy, the one on the NE 
being the far side. 

The most relevant observational facts are: a) the anoma- 
lous arms are made up of hot gas that has been shock ex- 
cited; b) the radio emission of the innermost region shows a 
clear jet (nearly in the N-S direction) , perpendicular to the 
maser disc around the black hole (The m aser disc is actu- 
ally ti pped down some 8° from edge on [see lHerrnstein et al.l 
l|l997l) ]. meaning that the northern jet is pointing toward us 
while the southern one is pointing away from us); c) the NW 
(SE) arm lies behind (in front of) the galaxy disc; and d) 
the inner section of the anomalous arms is roughly straight 
and it is roughly aligned with the major axis of the galaxy. 

From the several scenarios that have been proposed to 
explain the observational facts, the model that seems to be 
favour ed by the latest obser vations in X-rays ([Yang et al.l 
l[2007l) ) is the one proposed bv lWilson et~afl (l200lf )T in which 
a jet moving away from the disc induces strong shocks in the 
disc plane giving rise to the anomalous arms. 

We have made new observations of the inner region of 
NGC 4258 with an integral field spectrograph. These new 
observations will provide new kinematic and morphological 
information that will permit us to put the available models 
to the test. 



2 OBSERVATIONS, AND DATA REDUCTION 
AND ANALYSIS 

The data analysed in this paper were obtained on 2002 
March 16 at the Observatorio del Roque de los Muchachos 
on the island of La P alma with the INTEGRAL fibre system 
([Arribas et al1ll998f ) in combination with the WYFFOS fi- 
bre spectrograph (Bingham et alj Il994l ) at the William 
Herschel Telescope. Weather conditions during observations 
were fairly good, with an average seeing of 1.3 ". The data 
were obtained with INTEGRAL standard bundles SB3 and 
SB2, and the WYFFOS spectrograph was equipped with a 
1200 groove mm -1 grating centred on 6247 A. With these 
settings, the field of view is 33.6" x 29.4" and 16.0" x 12.3" 
respectively for the SB3 and SB2 bundles. The spectral res- 
olution is 4.8 A (R w 1300) for the SB3 bundle and 2.8 
A (R w 2200) for the SB2 bundle. The observed spectral 
range (5600-6850 A) contains a few strong interstellar emis- 
sion lines (Ha, [NII]AA6548,6584, [SII]AA6716, 6731), as well 
as several stellar absorption lines (NaD doublet and weaker, 
mostly blended Fel and Cal lines). Three exposures of 1200 s 
were taken for each fibre bundle. Standard data reduction for 
2D spectroscopic data (bias correction, flat-fielding, extrac- 
tion of spectra, wavel ength calibrat i on, sky subtraction, etc.) 
was p erformed (e.g. lArribas et alj l[l99ll ). iMediavilla et alj 
(1992)). The final wavelength calibration has uncertainties 
of 0.15 A or 7 km/s at 6500 A. 

In order to separate the stellar and interstellar contri- 
bution to the spectra, each spectrum is fitted using a syn- 
thetic stellar population as a template. The regions in the 
spectrum which are obviously contaminated by interstellar 
emission and/or absorption (namely Ha, [NIIJAA6548, 6584, 
[SII]AA6716,6731 and NaD) are not included in the fit. Af- 
ter some testing, we have finally used a population of so- 
lar metallicity with a n age of 1 Gy r from the synthetic 
population library by Delg ado et al. I l|2005h . which covers 
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Figure 1. Image of the inner region obtained with the F658N 
filter with the ACS on the HST. Overplotted contours show our 
observations with the same (synthetic) filter. 



the observed wavelength range at a very high spectral res- 
olution (0.3 A). Using populations with a slightly different 
age and/or metallicity do not seriously affect the results. 
We would like to point out that we do not intend a thor- 
ough analysis of the stellar populations in this region, but 
just to be able to subtract the stellar contribution from 
the observed spectra. As a by-product we obtain a stel- 
lar velocity, velocity dispersion and a global strength for 
the stellar lines at each fibre position. More details on this 
pro cedure and the general redu ctio n process can be found 
in ([Castillo-Morales et aT]|2007l ) and I Jimenez- Vicente et all 
(|2007h . 

Maps for several spectral features are built from the 
extracted spectra by mapping the extracted values at the 
fibre positions to a regular grid. The final images have 33 
x 29 pixels (with a size of 0.95 arcsec/pixel) and 45 x 34 
pixels (with a size of 0.35 arcsec/pixel) for the SB3 and 
SB2 bundles respectively. The resulting FOV are theferore 
of 31.35" x 27.55" for the SB3 and 15.75" x 11.9" for the SB2. 
Finally, astrometry was calculated for the maps by using the 
position of the nucleus and bright HII regions in the FOV. 



3 MORPHOLOGY AND KINEMATICS 

The optical images of this galaxy show very clearly that 
heavy absorption is taking place in t he (near side) SW part 
of the galaxy (see for example fig. 6 in lCecil et alj (|2000h ). It 
is very important to bear this fact in mind when interpreting 
the observations presented in this paper, particularly for the 
different effect it will have on the emission and absorption 
lines. 

The observed spectra show a few interesting emis- 
sion lines from the ISM, namely Ha, [Nil] AA6548, 6584, 
[SII]AA6716,6731 and [OI]AA6300, 6363. 

The distribution of Ha intensity from the SB3 and SB2 
bundles matches exceptionally well the features seen in the 
HST image taken by the ACS with the narrow F658N filter 
(see Fig. [!}. In order to compare both, we have generated 
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Figure 2. Continuum-free Ha image of the innermost region of 
NGC 4258 as observed with the SB2 fibre bundle of INTEGRAL. 
The circle on the upper right marks the defective fibre #167. 
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Figure 3. Star velocity field. The region in the SW is masked due 
to low S/N. Major axis is indicated. Contour levels are marked 
every 20 km/s. Systemic velocity of 450 km/s is marked as a 
thicker contour. Blue is approaching and red is receding. 



a synthetic F658N image from our data by integrating the 
spectra with the wavelengths weighted with the transmission 
curve of that filter. The image and contours show the above- 
mentioned strong absorption to the SW of the nucleus. In 
the (continuum-free) Ha image from the SB2 (with higher 
sp atial resolut i on), t he northern nuclear loop first detected 
bv lCecil et ail (2000) stands out even more clearly than in 
images from the HST (see Fig. [2]| . In our continuum-free Ha 
image, the contrast b etween the E and W sides of the loop 
seems weaker than in ICecil et al. | (|2000h . 

There is even a hint of a southern counterpart of this 
loop in the asymmetric morphology below the nucleus, al- 
though severe obscuration by dust in this region makes it 
impossible to confirm. The nature of the loop(s) remains 
unclear. A possible origin is that they are gas bubbles orig- 
inating from t he interaction of th e jet with the dense ISM, 
as proposed bv lCecil et al.l l|l992h . The velocity structure of 
the ionized gas, showing an asymmetry roughly along the 
N-S line, supports a connection with the jet. On the other 
hand, the line ratios at those locations are not indicative of 
shocks. Thus, the possibility of the lo op(s) being t he ed ges 
of an ionization cone, as proposed in ICecil et al.l (|2000l ) is 
also very attractive. 

The velocity field in Fig. [3] shows the rotation pattern 
of the stars. A big area in the SW part of the FOV has been 
masked because the low S/N in the spectra induced by the 
high extinction in this region does not allow us to calculate 
a reliable velocity. 

Even for the SB3 bundle, the FOV is quite small, and 
the velocity field of the ionized gas is too distorted to calcu- 
late a rotation model for the disc. We will therefore use the 
velocity field for the stars to subtract the rotating pattern 
of the galaxy from the gas components in order to gener- 
ate residual velocity fields. Obviously, there will be differ- 
ences between the rotating pattern of the stars and the gas 
(namely the asymmetric drift). We do expect these differ- 
ences to be small enough for our purposes and, in any case, 
we will be very cautious when interpreting the residual ve- 
locity fields. 



3.1 Neutral gas distribution and kinematics 

In this section we address the distribution and kinematics of 
the interstellar neutral gas. In Fig. [3] we present the distribu- 
tion of the neutral gas detected by means of the equivalent 
width of the NaD doublet. In Fig. [5] we show the residual 
velocities obtained by subtracting the stellar velocity map 
from the NaD velocity map. In the map of the equivalent 
width we can see a strip parallel to the major geometrical 
axis of the disc starting at the nucleus and extending some 
15 arcsec to each side of the nucleus. This strip nicely traces 
the neutral gas of the near side of the galaxy that is ab- 
sorbing the light coming from the stellar disc. Of course, 
this is not all the neutral gas, but only the one that can 
be detected in absorption (we cannot see the neutral gas 
on the far side of the galaxy, nor can we see the neutral 
gas in the more extinguished region on the nearest side for 
the same reason). Even with this obvious limitation, the de- 
tected distribution of cold neutral gas roughly resembles the 
distribution of molecular g as as seen in CO observations (cf. 
ISawada-Satoh et alj (120071 )'). 

Figure [5] shows the velocity of the neutral gas with 
respect to the underlying disc obtained by subtracting 
the stellar velocity map from the velocity map of the 
neutral gas. Again, the morphology of the map resem- 
bles the residual maps ob tained from CO observations by 
Sawada-S atoh et alj (|2007T > . but in our case the residual ve- 
locities are much higher. The detected neutral gas has huge 
residual velocities in the —120 to —3 70 km/s range with 
respec t to the stellar disc. Although ISawada-Satoh et alj 
(2007) have proposed a warped disc with expansion to ex- 
plain their data, we favour an alternative scenario to explain 
our detected morphology and kinematics of the neutral gas. 
In our opinion our data seem to indicate that the detected 
neutral gas (which is on the near side of the galaxy) it is 
most probably moving radially outwards from the galaxy 
centre. We will elaborate a bit more on this scenario later 
on. 
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Figure 4. Map of the equivalent width of the interstellar NaD. 
Regions with an EW lower than 0.8 Angstrom have been masked 
out. Brightest spot has an EW of 3.2. The major axis is indicated. 
The contours show synthetic image with F658N filter. 
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Figure 5. Map of the residual velocity field of the interstellar 
neutral gas. Velocity values are shown in km/s in the colour bar 
at the bottom. Contours show synthetic image with the F658 
filter. 



3.2 Kinematic anomalies and physical conditions 
of the ionized gas 

Comparing the velocity maps of the ionized gas (see Fig. [SJ 
with the stellar velocity field, we can see two conspicuous 
protrusions: one with receding velocities on the approaching 
side, and vice versa, roughly aligned at PA 145°. This strong 
disturbance of the velocity field is caused by a kinemati- 
cally distinct system of ionized gas that deviates from simple 
galactic rotation. We found direct evidence of this ionized 
gas of peculiar kinematics in the collection of spectra taken 
with the smaller SB2 bundle that provides higher spectral 



Figure 6. Map of the residual velocity field of the ionized gas. 
Velocity values are shown in km/s in the colour bar at the bot- 
tom. The thick continuous and dashed contours show regions with 
residuals higher than 50 km/s and lower than —50 km/s respec- 
tively. The major axis is indicated by the white line. 



resolution. We find many spectra in which the ionized lines 
split into two or more components with separations, in some 
cases, of mo re than 200 km/s. This line splitting was first 
reported by iRubin fc Graham! dl99Ch and an analysis was 
afterwards performed by Cecil et all 1 19921 ). who detected a 
braided jet structure. 

In order to attempt the identification of different 
gaseous systems in the area covered by our obse rvations, we 
have m ade a multi-component fit in the line of ICecil et al.l 
( 1992) based on the following assumptions: (i) there are two 
main kinematic systems of ionized gas, (ii) the broadening of 
the two components may be different, and (iii) the velocity 
of the two components is fairly different. Thus, we simul- 
taneously fit the two components to the N[II], SII and Ha 
emission lines. We accept the presence of two components 
when the goodness of the fit exceeds the goodness of the fit 
to a single component by 30%. Despite the intrinsic difficul- 
ties of this process and the peculiarities of this galaxy (e.g. in 
many cases the lines are not completely split, or there are ob- 
viously more than two components) which prevent the deter- 
mination of reliable velocity maps for both kinematic com - 
ponents (as was achieved by Ijimenez- Vicente et all (|2007T ) 
for a different object), it has been possible to identify con- 
sistently two main kinematic components at many locations 
throughout the FOV. The kinematics of the first one, which 
we identify with the gas in the galaxy disc, are not very dif- 
ferent from that of the stars. The other component, shows 
most peculiar, nearly counter-rotating kinematics. The lo- 
cations of the spectra where this anomalous component is 
identified are represented in Figure [7] We have chosen not 
to show the intensities, but only a sketch of the locations 
where two distinct kinematic components have been distin- 
guished. Intensities and velocities in both components may 
be seriously affected by the difficulties mentioned above, but 
locations where there is more than one component can be 
more reliably detected. It must be pointed out that there 
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Figure 7. Contour of the regions where a second kinematic com- 
ponent can be found in the ionized gas. The dashed lines show a 
sketch of the tentative bipolar structure. The major axis is indi- 
cated. 



still might be some locations where there are two or more 
components which have not been detected by our procedure, 
particularly in the more extinguished region to the south. 
This figure shows that this perturbed ionized gas component 
has a bipolar morphology. The northern region is somewhat 
inclined westwards, roughly between the (nearly vertical) 
axis defined by the jet and the major axis, and the southern 
part, is roughly eastwards from the (N-S) jet axis. The ex- 
istence of a southern counterpart for the northern Loop as 
mentioned in Section [3] is also hinted by this figure. 

In Figs. [9] and [10] the maps corresponding to the 
[NII]/Hq, [SII]/Hq and OI/Ha emission line ratios are pre- 
sented. According to these maps the central region has rel- 
atively low [NII]/Ha, [SII]/Hq and OI/Hq emission line ra- 
ti os. Taking into ac count the high [OIII]/H/3 measurements 
of lCecil et all (|l995h . our maps are consistent with photoion- 
ization as the main ionizing source for the central region, 
most probably induced by the Seyfert nucleus. Surrounding 
this central region towards the west there are several re- 
gions of very high [SII]/H« and OI/Ha emission line ratios 
forming an elongated strip almost parallel to the major ge- 
ometric axis. There is also a region with high ratios on the 
other side of the nucleus, suggesting a ring-like structure for 
the shocked gas. Shocks are the main cause of ionization 
of the gas in this strip/ring (front of shock). The border of 
the ionized gas component with peculiar kinematics overlaps 
with this shock front in the NW. A lthough shock-exci t ed gas 
ha s been prev i ously reported (cf. iRubin fe Graham! (| 19901 ) 
or ICecil et all jl995h ) we have been able to map the distri- 
bution of this shock ionized gas. 
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Figure 8. Map of the [NII]/Ha ratio with the SB3 bundle. Con- 
tours show synthetic image with the F658N filter. The major axis 
is indicated by the white line. 
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Figure 9. Map of the [SII]/Ha ratio. Contours show synthetic 
image with the F658N filter. The major axis is indicated by the 
white line. 



ISawada-Satoh et al.l (|2007f)); and (b) the bip olar hypershell 
model as described bv lSofue fc Vogleri l|200lf ). Both models 
may be supported by the present observations. 



4 DISCUSSION 

Among the variety of models proposed to describe the pe- 
culiar features of NGC 4258 that we have mentioned in 
Section [1] we will d iscuss two: (a) the model described by 
IWilson et alj |2001), which explains the anomalous arms 
as dense gas in the disc shocked by gas driven onto the 
disc by the out-of-plane cocoon generated by the jets, and 
which seems favoured by the latest X-ray observations (cf. 



4.1 Jet-driven shocks in the galaxy disc 

In Fig. [7] we can see that the bipolar component of the 
ionized gas with peculiar kinematics lies roughly in the re- 
gion between the N-S direction defined by the jet axis and 
the geometric major axis. This is the expected morphol- 
ogy in the scenario in which a cocoon of gas driven out of 
the plane by the jet is being blasted against the ISM disc 
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Figure 10. Map of the OI/Ha ratio. The contours show the 
synthetic image with the F658N filter. The major axis is indicated 
by the white line. 

(see I Wilson et all (|200lh and references therein). The bipo- 
lar component is mainly photoionized (most probably by 
the Seyfert nucleus) but its western border overlaps with 
the front shock detected in the [SII]/Hq (see Fig. [9]). This 
supports the notion that photoionized gas driven by the 
jet component is shocking with the ISM disc forming a hot 
cocoon/interface. According to this model, the front shock 
could trace the inner part of the anomalous arms. To com- 
plete the picture, the observed kinematics of the bipolar 
component of ionized gas (towards the blue in the NW and 
towards the red in the SE) will fit in this scheme if the co- 
coon/interface is back-flowing. Qualitatively, the presence 
of radial outflows in the neutral gas blasted by the cocoon 
would also support this model. 

4.2 The bipolar hypershell model 

ISofud (|l980h and ISofue fc Voglerl (|200ll ) interpreted the 
anomalous arms as the enhanced emission from the edges 
of a dumbbell-shaped shell that is symmetric with respect 
to the plane of NGC 4258. The presence of a wind creat- 
ing this structure is strongly supported by the detection 
of a very strong blueshift of the neutral gas emitters. This 
result is commonly accepted as solid evidence for a galac- 
tic wind, which is also supported by the existence of ion- 
ized gas systems of peculiar kinematics and the front shock. 
The main objection to Sofue's bip olar hypershell model 
(Ivan Albada fc van der HulstJ |l982l )) is the lack of symme- 
try of the hypershells. Although in Sofue's model the hyper- 
shells are generated by a starburst, the emission line ratios 
indicate that in NGC 4258 the entrai ning gas could be pho - 
toionized by the AGN. According to IVeilleux et alJ (|2005l ) 
when the source powering the galactic wind is an AGN, 
asymmetries in the expanding shell are expected. This is 
particularly true in the case of NGC 4258, where the out- 
of-plane jet would intersect the bipolar hypershell in the NS 
(projected) direction. The interaction between jet and hy- 



pershell would destroy that part of the hypershell bubbles 
(or would prevent it to form if the jet was in action when the 
shell was developing). Therefore, in this scenario, the anoma- 
lous arms would be the surviving parts of the hypershell. 
Indeed, the morphology of the rad io e mission as shown in 
Ivan Albada fc van der HulstJ l | 19821 ) and ICecil et al.l l|200Ch 
strongly suggests this interpretation. The sharp outer edges 
of the anomalous arms with the diffuse emission on the inside 
is naturally explained in this scenario. The observational fact 
that the SE arm is behind the ga l axy di sc while the NW arm 
is m front of it (cf. IWilson et all |200ll )') also fits very natu- 
rally into this model. Finally, the nearly straight morphology 
of the radio emission in the innermost region shows that the 
jet is clearly dominant in this region, in contrast with the 
outermost region, which is shell dominated. 



5 CONCLUSIONS 

We have performed two-dimensional spectroscopic observa- 
tions in the central region of NGC 4258 and have used these 
new data to test the available models for this peculiar ob- 
ject. The main results derived from the analysis of these new 
2D spectroscopy observations of NGC 4258 are: 

(i) The interstellar NaD absorption lines are strongly 
blueshifted. Comparison of the 2D velocity fields of the NaD 
with the stellar velocity field reveal that the blueshifted ve- 
locities are quite high, in the range —350 km/s to —120 
km/s, and occur everywhere in the detected neutral gas on 
the near side of the galaxy disc (that is, the one that we can 
observe in absorption). 

(ii) The velocity map of the ionized gas appears to be 
strongly distorted. There is a component of ionized gas of 
bipolar morphology with peculiar, almost counter-rotating, 
kinematics. 

(ill) The [SII]/Hq and OI/Ha emission line ratio maps 
show a thin ring-like region of high values indicative of 
shocks surrounding the galaxy nucleus. The border of the 
bipolar ionized gas component overlaps this front shock on 
the near side of the galaxy. 

The geometry of the bipolar ionized gas component and 
its partial overlapping with the front shock could fit into 
the model very well, in that hot gas driven by the jet blasts 
the ISM disc, giving rise to the anomalous arms. However, 
the radial motions of the ISM seem hard to explain by this 
model alone. This outflow, could be naturally interpre ted in 
terms of th e presence of the bipola r hypershell model (|Sofuel 
(|l980h and ISofue fc Voglerl (120011 )) in that the anomalous 
arms would be the enhanced emission from the edges of the 
shells. The lack of symmetry of the bubbles (lack of pres- 
ence of the anomalous arms in the NE and SW) would be 
explained by the disruption caused by the jet in the regions 
where it intersects the bubbles. 
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